Feedback cooling of a single trapped ion 
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Based on a real-time measurement of the motion of a single ion in a Paul trap, we demonstrate its 
electro-mechanical cooling below the Doppler limit by homodyne feedback control (cold damping). 
The feedback cooling results are well described by a model based on a quantum mechanical Master 
Equation. 
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Quantum optics, and more recently mesoscopic con- 
densed matter physics, have taken a leading role in re- 
alizing individual quantum systems, which can be mon- 
itored continuously in quantum limited measurements, 
and at the same time can be controlled by external fields 
on time scales fast in comparison with the system evolu- 
tion. Examples include cold trapped ions and atoms Ij, 
cavity QED 0- [J EH3 and nanomechanical systems g . 
This setting opens the possibility of manipulating indi- 
vidual quantum systems by feedback, a problem which is 
not only of a fundamental interest in quantum mechan- 
ics, but also promises a new route to generating interest- 
ing quantum states in the laboratory. First experimental 
efforts to realize quantum feedback have been reported 
only recently. While not all of them may qualify as quan- 
tum feedback in a strict sense,feedback has been applied 
to various quantum systems 0, 0, la M> LLil 111! ■ On the 
theory side, this has motivated during the last decade the 
development of a quantum feedback theory |12U13| ]. where 
the basic ingredients are the interplay between quantum 
dynamics and the back-action of the measurement on 
the system evolution. In this letter we report a first 
experiment to demonstrate quantum feedback control, 
i.e. quantum feedback cooling, of a single trapped ion 
by monitoring the fluorescence of the laser driven ion in 
front of a mirror. We establish a continuous measurement 
of the position of the ion which allows us to act back in 
a feedback loop demonstrating "cold damping" [Til fl5| . 
We will show that quantum control theory based on a 
quantum optical modelling of the system dynamics and 
continuous measurement theory of photodetection pro- 
vides a quantitative understanding of the experimental 
results. 

We study a single 138 Ba + ion in a miniature Paul trap 
which is continuously laser-excited and laser-cooled to 
the Doppler limit on its S 1( / 2 to P1/2 transition at 493 nm, 
as outlined in Fig. 1. The ion is driven by a laser near 
the atomic resonance, and the scattered light is emitted 
both into the radiation modes reflected by the mirror, as 
well as the other (background) modes of the quantized 



light field [16j . Light scattered into the mirror modes can 
either reach the photodetector directly, or after reflection 
from the mirror. From the resulting interference the mo- 
tion of the ion (its projection onto the ion-mirror axis) 
is detected as a vibrational sideband in the fluctuation 
spectrum of the photon counting signal Of the three 
sidebands at about (1,1.2,2.3) MHz, corresponding to the 
three axes of vibration, we observe the one at v = 1 MHz. 
It has a width T « 400 Hz and is superimposed on the 
background shot noise generated by the photon counting 
process. 

Our goal is to continuously read the position of the ion, 
and feed back a damping force proportional to the mo- 
mentum to achieve feedback cooling. For a weakly driven 
atom the emitted light is dominantly elastic scattering at 
the laser frequency, and the information on the motion of 
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FIG. 1: A single 138 Ba + ion in a Paul trap (parabola) is laser- 
excited and -cooled on its S 1 / 2 to P1/2 transition at 493 nm. 
A retro-reflecting mirror 25 cm away from the trap and a 
lens (not shown) focus back the fluorescence onto the ion. 
The resulting interference fringes with up to 73% contrast are 
observed by a photomultiplier (PMT 1). The ion's oscillation 
in the trap creates an intensity modulation of the PMT signal 
which is observed as a sideband on a spectrum analyser (rfsa) 
[r?| . For feedback cooling, the sideband signal is filtered, 
phase-shifted, and applied to the ion as a voltage on a trap 
electrode. 
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the ion is encoded in the sidebands of the scattered light, 
displaced by the trap frequency v. For an ion trapped 
in the Lamb-Dicke regime (rjy/N <C 1, with the Lamb 
Dicke parameter rj = 27ra /A ~ 0.07, ao the r.m.s. size 
of the trap ground state, N the mean excitation num- 
ber of the motional oscillator) the motional sidebands 
are suppressed by rj relative to the elastic component at 
the laser frequency. Thus the light reaching the detector 
will be the sum of the elastic component and weak side- 
bands |l7j . a situation reminiscent of homodyne detec- 
tion, where a strong oscillator beats with the signal field 
to provide a homodyne current at the detector. This 
physical picture allows us to formulate the continuous 
readout of the position of the ion as well as the quan- 
tum feedback cooling in the well-developed language of 
homodyne detection and quantum feedback. 

The homodyne current at the photodetector (see 
Fig. 1) with the (large) signal from the elastic light scat- 
tering subtracted has the form 

Io{t)=ir ) {z) c {t) + s jl^t) . (1) 

The first term is proportional to the conditioned expecta- 
tion value of the position of the trapped ion, (z) c (t), and 
the second term is a shot noise contribution with Gaus- 
sian white noise We have defined z = a + a) = z/ao 
with a (of) destruction (creation) operator of the har- 
monic oscillator, and we have assumed that the trap cen- 
ter is placed at a distance L = n\/2+X/8 (n integer) from 
the mirror, corresponding to a point of maximum slope of 
the standing wave intensity of the mirror mode. The cur- 
rent I c (t) scales with 70c e, which is the light scattering 
rate into the solid angle (47re) of the mirror mode induced 
by the laser. The expectation value (-) c = Tr{- p c (t)} is 
defined with respect to a conditioned density operator 
p c (t), which reflects our knowledge of the motional state 
of the ion for the given history of the photocurrent. Ac- 
cording to the theory of homodyne detection, p c (t) obeys 
the Ito stochastic differential equation 0] 

dp c (t) — — iis{a) a, p c {t)]dt + £op c (t)dt 

+ ^ Hp c (t)dW(t) , (2) 

where H Pc (t) = (zp c (t) + p c (t)z - 2 (z) c (t)p c (t)) . The 
first line determines the unobserved evolution of the ion, 
including the harmonic motion in the trap with frequency 
v and the dissipative dynamics, Co, due to photon scat- 
tering. The latter is given by 

C p = T(N + l)V[a]p + TNV[^]p , (3) 

where we have defined the superoperator T)[c]p= epe' — 
(c^cp + pc { c)/2. The laser cooling rate T and the steady 
state occupation number, N = (aJa), can either be es- 
timated from the motional sidebands or deduced from 
the cooling laser parameters 0. In the present exper- 
iment, N sa 17 corresponds to the Doppler limit. The 



last term of Eq. J5J is proportional to the Wiener incre- 
ment dW(t) = £(t)dt and corresponds to an update of 
the observers knowledge about the system according to 
a certain measurement result I c (t). In summary, Eq. Q 
demonstrates that observation of the sidebands of the 
light scattered into the mirror mode provides us with in- 
formation of the position of the ion, while the system 
density matrix is updated according to Eq. (J2J). This is 
the basis for describing feedback control of the ion, as 
shown in the following. 

For feedback cooling, the vibrational sideband is ex- 
tracted with a bandpass filter of bandwidth B = 30 kHz, 
shifted by (— 7r/2), amplified, and the resulting output 
voltage is applied to an electrode which is close to the 
trap inside the vacuum. Thereby we create a driving 
force which is proportional to (and opposed to) the in- 
stantaneous velocity of the ion and which thus adds to 
the damping of its vibration. The overall gain of the 
feedback loop depends on many factors such as the fringe 
contrast of the interference, the PMT characteristics, etc. 
It is varied electronically by setting the amplification G of 
the final amplifier in the loop. We can also set the phase 
to other values than the optimum of (— 7r/2), which will 
be used in order to compare experiment and theory. 

To analyze the result of the feedback we look at the 
changes in the sideband spectrum. The modified spectra 
require careful interpretation. The spectrum observed 
inside the feedback loop ("in-loop" or PMT1 in Fig. HJ 
shows not only the motion of the ion, but the sum of 
the motion and the shot noise. As the feedback cor- 
relates these fluctuations, a reduction of the signal be- 
low the shot noise level may occur, similar in appear- 
ance to signatures of squeezed light. This effect is known 
as "squashing" or "anti-correlated state of light" 

in an opto-electronic loop [2(|, and it does not consti- 
tute a quantum mechanical squeezing of the fluctuations 
[2ll E3 | . The effect on the motion can only be reliably 
detected by splitting the optical signal before it is mea- 
sured and recording it outside the feedback loop with a 
second PMT ("out-loop" in Fig.^l, whose shot noise is 
not correlated with the motion. 

In Fig. |21 we show spectra recorded with the spectrum 
analyzer, and measured outside and inside the feedback 
loop. The first curve of each row, showing the largest 
sideband, is the one without feedback (gain G = 0). The 
other two curves are recorded with the loop closed (gain 
values G — 1.3 and 8.7). The sub shot noise fluctuations 
inside the loop, when the ion is driven to move in anti- 
phase with the shot noise, are clearly visible in Fig. Elf). 
The main cooling results are curves (b) and (c) which 
show the motional sideband reduced in size and broad- 
ened, indicating a reduced energy (proportional to the 
area under the curve) and a higher damping rate (the 
width). From case (b) to (c) the area increases, as the 
injected and amplified shot noise overcompensates the 
increased damping. As shown below, in our model the 
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FIG. 2: Feedback cooling spectra. The vibrational sideband 
around v — 1 MHz is shown on top of the spectrally flat shot 
noise background, to which the spectra are normalized. The 
upper curves (a)-(c) are measured outside the feedback loop, 
while the lower curves (d)-(f) are the in-loop results. Spectra 
(a) and (d) are for laser cooling only, the other curves are 
recorded with feedback at the indicated gain values. The 
feedback phase is set to (— 7r/2). The gain values indicate the 
settings of the final amplifier in the feedback loop. 



incorporation of quantum feedback competing with laser 
cooling predicts such behavior, i.e., the existence of an 
optimal gain for maximal cooling (for a detailed descrip- 
tion cf. Hi). 

We model the effect of the feedback force acting on the 
ion by extending Eq. J5J with the feedback contribution, 



where (p) c (t) = Tr{p[i c (t)}, and p = i(a t — a) is the mo- 
mentum operator conjugated to z. The first term in this 
current therefore provides damping for the motion of the 
ion. The second term of Eq. (JSJ) describes the shot noise 
which passes through the electronic circuit and is fed back 
to the ion. The stochastic variable S(i) is Gaussian white 
noise on a timescale given by the inverse bandwidth B , 
whereby B ^> T implies that it is spectrally flat on the 
frequency range of the cooling dynamics. 

For a full record of the photocurrent I c (t), Eqs. (0 and 
(@J determine the evolution of the ion's motional state in 
presence of feedback. As it is impractical to keep track 
of the whole photocurrent in the experiment, we derive a 
master equation for the density operator averaged over all 
possible realisation of I c (t), n(t). Along the ways of the 
Wiseman-Milburn theory of quantum feedback f° r 
a phase shift of (— tt/2), we obtain the quantum feedback 
master equation [2^| 
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(6) 



The second and third term are the additional contribu- 
tions due to the feedback. The part linear in G induces 
damping of the motion of the ion, and the term quadratic 
in G describes the effect of the fed back noise leading to 
diffusion of the momentum. The competition between 
laser cooling, damping and the injected noise leads to 
the characteristic behavior of the steady state number 
expectation value 



[d Pc (t)]ib = -iGI !hiC (t - T )[z, p c (t)], 



(4) 



where Ifb, c {t) denotes the measured current after the 
feedback circuit. The latter is proportional to the voltage 
applied to the trap electrodes. All conversion factors be- 
tween the feedback current and the actual force applied 
on the ion are included in the overall gain G. The time 
delay r in the feedback loop preserves causality and is 
small compared to the fastest timescale z/ _1 of the mo- 
tion of the ion which allows us to consider the Markovian 
limit (t-> 0+). 

To obtain an expression for the feedback current 
I{b,c(t), we change into a frame rotating with the trap 
frequency v and define the density operator, p c (t) = 
exp(iisa^at)p c (t) exp(— iva^at), evolving on the (slow) 
cooling timescales. This is convenient due to the large 
separation between the timescale of the harmonic oscil- 



lations, 



and the timescale of laser and feedback cool- 



ing in our experiment. For our experimental parameters, 
v ^> B r, the feedback current for a phase shift of 
(-vr/2) has the form 



cos(j4), (5) 



N + r]-/G(2N - l)/2r + 7 G 2 /8r 
1 + 2777G/T ' 



(7) 



For small gain, damping dominates, and the energy of 
the ion is decreased below the Doppler limit. For higher 
gain, the diffusive term describing the noise fed back into 
the system overcompensates cooling, i.e., heats the ion. 
Consequently, for (— tt/2) feedback phase and optimal 
gain conditions the steady state energy is minimized. On 
the contrary, for a tt phase shift (— z) replaces p in the 
second term of Eq. the feedback force then merely 
induces a frequency shift, Au> = G777/2, but no damp- 
ing. Increasing the gain then always enhances the steady 
state number expectation value, i.e., the mean ion en- 
ergy. We now compare the theoretical predictions and 
the experimental results for the feedback phases of phys- 
ical interest, namely (— 7r/2) and 7r. 

The measured ion energy as a function of the feed- 
back electronic gain is shown in Fig. [21 On the first 
side, as expected, for a (—tt/2) feedback phase, cooling 
by more than 30% below the Doppler limit is achieved, 
while further increase of the gain drives the shot noise 
and therefore heats the motion of the ion. On the other 
side, a tt phase shift in the feedback loop does not yield 
any damping, in such conditions the motion of the ion is 
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FIG. 3: Steady-state energy of the cooled oscillator: measured 
sideband area, normalized to the value without feedback, ver- 
sus gain of the feedback loop, for (— 7r/2) (a) and n (b) feed- 
back phase. The curves are the model calculations. The gain 
axis is scaled to the experimental values of the electronic gain. 



only driven. This results in an increase of the measured 
sideband area (as shown in the insert of Figure EJ) , as 
well as a shift of the sideband center frequency (graph 
not shown) . Both cases demonstrate good agreement be- 
tween experiment and theory. Finally, let us stress that 
the optimal cooling rate is governed by the collection ef- 
ficiency of the fluorescence going into the mirror mode, 
e. In the experiments presented above, e ~ 1% leads 
to a decrease of the steady state occupation number iV 
from 17 to 12, while A « 3 can be reached for e « 15%, 
experimentally achievable for optimal optical coupling. 

To summarize, we have demonstrated real-time feed- 
back cooling of the motion of a single trapped ion. 
Electro-mechanical back-action based on a sensitive real- 
time measurement of the motion of the ion in the trap 
allowed us to cool one motional degree of freedom by 30% 
below the Doppler limit. Unlike with laser cooling, the 
presented method allows to cool one of the ion's motional 
mode without heating the two others, and our procedure 
can easily be extended to cooling all motional modes. 
The cooling process is shot-noise limited and the frac- 
tion of scattered photons recorded to observe the motion 
of the ion limits the ultimate cooling at optimum gain. 
The latter can yield a steady state occupation number 
N = 3 for realistic experimental conditions. In conclu- 
sion, our feedback scheme offers a possible way to very 
efficiently cool the motion of ions unsuitable for sideband 
cooling. 
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